The dynamics of the flavin, NADPH, and active site loops determine the mechanism of activation of class B flavin-dependent monooxygenases.

Flavin-dependent monooxygenases (FMOs) constitute a diverse enzyme family that catalyzes crucial hydroxylation, epoxidation, and Baeyer-Villiger reactions across various metabolic pathways in all domains of life. Due to the intricate nature of this enzyme family's mechanisms, some aspects of their functioning remain unknown. Here, we present the results of molecular dynamics computations, supplemented by a bioinformatics analysis, that clarify the early stages of their catalytic cycle. We have elucidated the intricate binding mechanism of NADPH and L-Orn to a class B monooxygenase, the ornithine hydroxylase from Aspergillus $$ Aspergillus $$ fumigatus $$ fumigatus $$ known as SidA. Our investigation involved a comprehensive characterization of the conformational changes associated with the FAD (Flavin Adenine Dinucleotide) cofactor, transitioning from the out to the in position. Furthermore, we explored the rotational dynamics of the nicotinamide ring of NADPH, shedding light on its role in facilitating FAD reduction, supported by experimental evidence. Finally, we also analyzed the extent of conservation of two Tyr-loops that play critical roles in the process.

Structural insights into the semiquinone form of human Cytochrome P450 reductase by DEER distance measurements between a native flavin and a spin labelled non-canonical amino acid.

The flavoprotein Cytochrome P450 reductase (CPR) is the unique electron pathway from NADPH to Cytochrome P450 (CYPs). The conformational dynamics of human CPR in solution, which involves transitions from a "locked/closed" to an "unlocked/open" state, is crucial for electron transfer. To date, however, the factors guiding these changes remain unknown. By Site-Directed Spin Labelling coupled to Electron Paramagnetic Resonance spectroscopy, we have incorporated a non-canonical amino acid onto the flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) domains of soluble human CPR, and labelled it with a specific nitroxide spin probe. Taking advantage of the endogenous FMN cofactor, we successfully measured for the first time, the distance distribution by DEER between the semiquinone state FMNH• and the nitroxide. The DEER data revealed a salt concentration-dependent distance distribution, evidence of an "open" CPR conformation at high salt concentrations exceeding previous reports. We also conducted molecular dynamics simulations which unveiled a diverse ensemble of conformations for the "open" semiquinone state of the CPR at high salt concentration. This study unravels the conformational landscape of the one electron reduced state of CPR, which had never been studied before.

The enzymes of the oxidative phase of the pentose phosphate pathway as targets of reactive species: consequences for NADPH production.

The pentose phosphate pathway (PPP) is a key metabolic pathway. The oxidative phase of this process involves three reactions catalyzed by glucose-6-phosphate dehydrogenase (G6PDH), 6-phosphogluconolactonase (6PGL) and 6-phosphogluconate dehydrogenase (6PGDH) enzymes. The first and third steps (catalyzed by G6PDH and 6PGDH, respectively) are responsible for generating reduced nicotinamide adenine dinucleotide phosphate (NAPDH), a key cofactor for maintaining the reducing power of cells and detoxification of both endogenous and exogenous oxidants and electrophiles. Despite the importance of these enzymes, little attention has been paid to the fact that these proteins are targets of oxidants. In response to oxidative stimuli metabolic pathways are modulated, with the PPP often up-regulated in order to enhance or maintain the reductive capacity of cells. Under such circumstances, oxidation and inactivation of the PPP enzymes could be detrimental. Damage to the PPP enzymes may result in a downward spiral, as depending on the extent and sites of modification, these alterations may result in a loss of enzymatic activity and therefore increased oxidative damage due to NADPH depletion. In recent years, it has become evident that the three enzymes of the oxidative phase of the PPP have different susceptibilities to inactivation on exposure to different oxidants. In this review, we discuss existing knowledge on the role that these enzymes play in the metabolism of cells, and their susceptibility to oxidation and inactivation with special emphasis on NADPH production. Perspectives on achieving a better understanding of the molecular basis of the oxidation these enzymes within cellular environments are given.

From a binding module to essential catalytic activity: how nature stumbled on a good thing.

Enzymes are complex macromolecules capable of catalyzing a wide variety of chemical reactions with high efficiency. Nonetheless, biological catalysis can be rudimentary. Here, we describe an enzyme that is built from a simple protein fold. This short protein sequence - almost a peptide - belongs to the ancient SH3 family of binding modules. Surprisingly, this binding module catalyzes the specific reduction of dihydrofolate using NADPH as a reducing cofactor, making this a dihydrofolate reductase. Too small to provide all the required binding and catalytic machinery on its own, it homotetramerizes, thus creating a large, central active site environment. Remarkably, none of the active site residues is essential to the catalytic function. Instead, backbone interactions juxtapose the reducing cofactor proximal to the target imine of the folate substrate, and a specific motion of the substrate promotes formation of the transition state. In this feature article, we describe the features that make this small protein a functional enzyme capable of catalyzing a metabolically essential reaction, highlighting the characteristics that make it a model for the evolution of primitive enzymes from binding modules.

Structure of the imine reductase from Ajellomyces dermatitidis in three crystal forms.

The NADPH-dependent imine reductase from Ajellomyces dermatitidis (AdRedAm) catalyzes the reductive amination of certain ketones with amine donors supplied in an equimolar ratio. The structure of AdRedAm has been determined in three forms. The first form, which belongs to space group P3121 and was refined to 2.01 Šresolution, features two molecules (one dimer) in the asymmetric unit in complex with the redox-inactive cofactor NADPH4. The second form, which belongs to space group C21 and was refined to 1.73 Šresolution, has nine molecules (four and a half dimers) in the asymmetric unit, each complexed with NADP+. The third form, which belongs to space group P3121 and was refined to 1.52 Šresolution, has one molecule (one half-dimer) in the asymmetric unit. This structure was again complexed with NADP+ and also with the substrate 2,2-difluoroacetophenone. The different data sets permit the analysis of AdRedAm in different conformational states and also reveal the molecular basis of stereoselectivity in the transformation of fluorinated acetophenone substrates by the enzyme.

The unusual chemical sequences of mammalian dihydropyrimidine dehydrogenase revealed by transient-state analysis.

Dihydropyrimidine dehydrogenase (DPD) catalyzes the reduction of the 5,6-vinylic bond of uracil and thymine with electrons from NADPH. The complexity of the enzyme belies the simplicity of the reaction catalyzed. To accomplish this chemistry DPD has two active sites that are ∼60Šapart, both of which house flavin cofactors, FAD and FMN. The FAD site interacts with NADPH, while the FMN site with pyrimidines. The distance between the flavins is spanned by four Fe4S4 centers. Though DPD has been studied for nearly 50years, it is only recently that the novel apects of its mechanism have been described. The primary reason for this is that the chemistry of DPD is not portrayed adequately by known descriptive steady-state mechanism categories. The highly chromophoric nature of the enzyme has recently been exploited in transient-state to document unexpected reaction sequences. Specifically, DPD undergoes reductive activation prior to catalytic turnover. Two electrons are taken up from NADPH and transmitted via the FAD and Fe4S4 centers to form the FAD•4(Fe4S4)•FMNH2 form of the enzyme. This form of the enzyme will only reduce pyrimidine substrates in the presence NADPH, establishing that hydride transfer to the pyrimidine precedes reductive reactivation that reinstates the active form of the enzyme. DPD is therefore the first flavoprotein dehydrogenase known to complete the oxidative half-reaction prior to the reductive half-reaction. Here we describe the methods and deduction that led to this mechanistic assignment.

Sequence and structure-guided discovery of a novel NADH-dependent 7ß-hydroxysteroid dehydrogenase for efficient biosynthesis of ursodeoxycholic acid.

7ß-Hydroxysteroid dehydrogenases (7ß-HSDHs) have attracted increasing attention due to their crucial roles in the biosynthesis of ursodeoxycholic acid (UDCA). However, most published 7ß-HSDHs are strictly NADPH-dependent oxidoreductases with poor activity and low productivity. Compared with NADPH, NADH is more stable and cheaper, making it the more popular cofactor for industrial applications of dehydrogenases. Herein, by using a sequence and structure-guided genome mining approach based on the structural information of conserved cofactor-binding motifs, we uncovered a novel NADH-dependent 7ß-HSDH (Cle7ß-HSDH). The Cle7ß-HSDH was overexpressed, purified, and characterized. It exhibited high specific activity (9.6 U/mg), good pH stability and thermostability, significant methanol tolerance, and showed excellent catalytic efficiencies (kcat/Km) towards 7-oxo-lithocholic acid (7-oxo-LCA) and NADH (70.8 mM-1s-1 and 31.8 mM-1s-1, respectively). Molecular docking and mutational analyses revealed that Asp42 could play a considerable role in NADH binding and recognition. Coupling with a glucose dehydrogenase for NADH regeneration, up to 20 mM 7-oxo-LCA could be completely transformed to UDCA within 90 min by Cle7ß-HSDH. This study provides an efficient approach for mining promising enzymes from genomic databases for cost-effective biotechnological applications.

Crystal structure of Leishmania donovani glucose 6-phosphate dehydrogenase reveals a unique N-terminal domain.

Since unicellular parasites highly depend on NADPH as a source for reducing equivalents, the pentose phosphate pathway, especially the first and rate-limiting NADPH-producing enzyme glucose 6-phosphate dehydrogenase (G6PD), is considered an excellent antitrypanosomatid drug target. Here we present the crystal structure of Leishmania donovani G6PD (LdG6PD) elucidating the unique N-terminal domain of Kinetoplastida G6PDs. Our investigations on the function of the N-domain suggest its involvement in the formation of a tetramer that is completely different from related Trypanosoma G6PDs. Structural and functional investigations further provide interesting insights into the binding mode of LdG6PD, following an ordered mechanism, which is confirmed by a G6P-induced domain shift and rotation of the helical N-domain. Taken together, these insights into LdG6PD contribute to the understanding of G6PDs' molecular mechanisms and provide an excellent basis for further drug discovery approaches.

Unravelling the regulation pathway of photosynthetic AB-GAPDH.

Oxygenic phototrophs perform carbon fixation through the Calvin-Benson cycle. Different mechanisms adjust the cycle and the light-harvesting reactions to rapid environmental changes. Photosynthetic glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is a key enzyme in the cycle. In land plants, different photosynthetic GAPDHs exist: the most abundant isoform is formed by A2B2 heterotetramers and the least abundant by A4 homotetramers. Regardless of the subunit composition, GAPDH is the major consumer of photosynthetic NADPH and its activity is strictly regulated. While A4-GAPDH is regulated by CP12, AB-GAPDH is autonomously regulated through the C-terminal extension (CTE) of its B subunits. Reversible inhibition of AB-GAPDH occurs via the oxidation of a cysteine pair located in the CTE and the substitution of NADP(H) with NAD(H) in the cofactor-binding site. These combined conditions lead to a change in the oligomerization state and enzyme inhibition. SEC-SAXS and single-particle cryo-EM analysis were applied to reveal the structural basis of this regulatory mechanism. Both approaches revealed that spinach (A2B2)n-GAPDH oligomers with n = 1, 2, 4 and 5 co-exist in a dynamic system. B subunits mediate the contacts between adjacent tetramers in A4B4 and A8B8 oligomers. The CTE of each B subunit penetrates into the active site of a B subunit of the adjacent tetramer, which in turn moves its CTE in the opposite direction, effectively preventing the binding of the substrate 1,3-bisphosphoglycerate in the B subunits. The whole mechanism is made possible, and eventually controlled, by pyridine nucleotides. In fact, NAD(H), by removing NADP(H) from A subunits, allows the entrance of the CTE into the active site of the B subunit, hence stabilizing inhibited oligomers.

Allosteric role of a structural NADP+ molecule in glucose-6-phosphate dehydrogenase activity.

Human glucose-6-phosphate dehydrogenase (G6PD) is the main cellular source of NADPH, and thus plays a key role in maintaining reduced glutathione to protect cells from oxidative stress disorders such as hemolytic anemia. G6PD is a multimeric enzyme that uses the cofactors ß-D-glucose 6-phosphate (G6P) and "catalytic" NADP+ (NADP+c), as well as a "structural" NADP+ (NADP+s) located ∼25 Å from the active site, to generate NADPH. While X-ray crystallographic and biochemical studies have revealed a role for NADP+s in maintaining the catalytic activity by stabilizing the multimeric G6PD conformation, other potential roles for NADP+s have not been evaluated. Here, we determined the high resolution cryo-electron microscopy structures of human wild-type G6PD in the absence of bound ligands and a catalytic G6PD-D200N mutant bound to NADP+c and NADP+s in the absence or presence of G6P. A comparison of these structures, together with previously reported structures, reveals that the unliganded human G6PD forms a mixture of dimers and tetramers with similar overall folds, and binding of NADP+s induces a structural ordering of a C-terminal extension region and allosterically regulates G6P binding and catalysis. These studies have implications for understanding G6PD deficiencies and for therapy of G6PD-mediated disorders.

Conformational changes in the catalytic region are responsible for heat-induced activation of hyperthermophilic homoserine dehydrogenase.

When overexpressed as an immature enzyme in the mesophilic bacterium Escherichia coli, recombinant homoserine dehydrogenase from the hyperthermophilic archaeon Sulfurisphaera tokodaii (StHSD) was markedly activated by heat treatment. Both the apo- and holo-forms of the immature enzyme were successively crystallized, and the two structures were determined. Comparison among the structures of the immature enzyme and previously reported structures of mature enzymes revealed that a conformational change in a flexible part (residues 160-190) of the enzyme, which encloses substrates within the substrate-binding pocket, is smaller in the immature enzyme. The immature enzyme, but not the mature enzyme, formed a complex that included NADP+, despite its absence during crystallization. This indicates that the opening to the substrate-binding pocket in the immature enzyme is not sufficient for substrate-binding, efficient catalytic turnover or release of NADP+. Thus, specific conformational changes within the catalytic region appear to be responsible for heat-induced activation.

Oleic acid based experimental evolution of Bacillus megaterium yielding an enhanced P450 BM3 variant.

BACKGROUND: Unlike most other P450 cytochrome monooxygenases, CYP102A1 from Bacillus megaterium (BM3) is both soluble and fused to its redox partner forming a single polypeptide chain. Like other monooxygenases, it can catalyze the insertion of oxygen unto the carbon-hydrogen bond which can result in a wide variety of commercially relevant products for pharmaceutical and fine chemical industries. However, the instability of the enzyme holds back the implementation of a BM3-based biocatalytic industrial processes due to the important enzyme cost it would prompt. RESULTS: In this work, we sought to enhance BM3's total specific product output by using experimental evolution, an approach not yet reported to improve this enzyme. By exploiting B. megaterium's own oleic acid metabolism, we pressed the evolution of a new variant of BM3, harbouring 34 new amino acid substitutions. The resulting variant, dubbed DE, increased the conversion of the substrate 10-pNCA to its product p-nitrophenolate 1.23 and 1.76-fold when using respectively NADPH or NADH as a cofactor, compared to wild type BM3. CONCLUSIONS: This new DE variant, showed increased organic cosolvent tolerance, increased product output and increased versatility in the use of either nicotinamide cofactors NADPH and NADH. Experimental evolution can be used to evolve or to create libraries of evolved BM3 variants with increased productivity and cosolvent tolerance. Such libraries could in turn be used in bioinformatics to further evolve BM3 more precisely. The experimental evolution results also supports the hypothesis which surmises that one of the roles of BM3 in Bacillus megaterium is to protect it from exogenous unsaturated fatty acids by breaking them down.

Investigating effect of mutation on structure and function of G6PD enzyme: a comparative molecular dynamics simulation study.

Several natural mutants of the human G6PD enzyme exist and have been reported. Because the enzymatic activities of many mutants are different from that of the wildtype, the genetic polymorphism of G6PD plays an important role in the synthesis of nucleic acids via ribulose-5-phosphate and formation of reduced NADP in response to oxidative stress. G6PD mutations leading to its deficiency result in the neonatal jaundice and acute hemolytic anemia in human. Herein, we demonstrate the molecular dynamics simulations of the wildtype G6PD and its three mutants to monitor the effect of mutations on dynamics and stability of the protein. These mutants are Chatham (A335T), Nashville (R393H), Alhambra (V394L), among which R393H and V394L lie closer to binding site of structural NADP+. MD analysis including RMSD, RMSF and protein secondary structure revealed that decrease in the stability of mutants is key factor for loss of their activity. The results demonstrated that mutations in the G6PD sequence resulted in altered structural stability and hence functional changes in enzymes. Also, the binding site, of structural NADP+, which is far away from the catalytic site plays an important role in protein stability and folding. Mutation at this site causes changes in structural stability and hence functional deviations in enzyme structure reflecting the importance of structural NADP+ binding site. The calculation of binding free energy by post processing end state method of Molecular Mechanics Poisson Boltzmann SurfaceArea (MM-PBSA) has inferred that ligand binding in wildtype is favorable as compared to mutants which represent destabilised protein structure due to mutation that in turn may hinder the normal physiological function. Exploring individual components of free energy revealed that the van der Waals energy component representing non-polar/hydrophobic energy contribution act as a dominant factor in case of ligand binding. Our study also provides an insight in identifying the key inhibitory site in G6PD and its mutants which can be exploited to use them as a target for developing new inhibitors in rational drug design.

Enhancing cofactor regeneration of cyanobacteria for the light-powered synthesis of chiral alcohols.

Cyanobacteria Synechocystis sp. PCC 6803 was exploited as green cell factory for light-powered asymmetric synthesis of aromatic chiral alcohols. The effect of temperature, light, substrate and cell concentration on substrate conversions were investigated. Under the optimal condition, a series of chiral alcohols were synthesized with conversions up to 95% and enantiomer excess (ee) > 99%. We found that the addition of Na2S2O3 and Angeli's Salt increased the NADPH content by 20% and 25%, respectively. As a result, the time to reach 95% substrate conversion was shortened by 12 h, which demonstrated that the NADPH regeneration and hence the reaction rates can be regulated in cyanobacteria. This blue-green algae based biocatalysis showed its potential for chiral compounds production in future.

Rossmann-toolbox: a deep learning-based protocol for the prediction and design of cofactor specificity in Rossmann fold proteins.

The Rossmann fold enzymes are involved in essential biochemical pathways such as nucleotide and amino acid metabolism. Their functioning relies on interaction with cofactors, small nucleoside-based compounds specifically recognized by a conserved ßαß motif shared by all Rossmann fold proteins. While Rossmann methyltransferases recognize only a single cofactor type, the S-adenosylmethionine, the oxidoreductases, depending on the family, bind nicotinamide (nicotinamide adenine dinucleotide, nicotinamide adenine dinucleotide phosphate) or flavin-based (flavin adenine dinucleotide) cofactors. In this study, we showed that despite its short length, the ßαß motif unambiguously defines the specificity towards the cofactor. Following this observation, we trained two complementary deep learning models for the prediction of the cofactor specificity based on the sequence and structural features of the ßαß motif. A benchmark on two independent test sets, one containing ßαß motifs bearing no resemblance to those of the training set, and the other comprising 38 experimentally confirmed cases of rational design of the cofactor specificity, revealed the nearly perfect performance of the two methods. The Rossmann-toolbox protocols can be accessed via the webserver at https://lbs.cent.uw.edu.pl/rossmann-toolbox and are available as a Python package at https://github.com/labstructbioinf/rossmann-toolbox.

Reductive inactivation of the hemiaminal pharmacophore for resistance against tetrahydroisoquinoline antibiotics.

Antibiotic resistance is becoming one of the major crises, among which hydrolysis reaction is widely employed by bacteria to destroy the reactive pharmacophore. Correspondingly, antibiotic producer has canonically co-evolved this approach with the biosynthetic capability for self-resistance. Here we discover a self-defense strategy featuring with reductive inactivation of hemiaminal pharmacophore by short-chain dehydrogenases/reductases (SDRs) NapW and homW, which are integrated with the naphthyridinomycin biosynthetic pathway. We determine the crystal structure of NapW·NADPH complex and propose a catalytic mechanism by molecular dynamics simulation analysis. Additionally, a similar detoxification strategy is identified in the biosynthesis of saframycin A, another member of tetrahydroisoquinoline (THIQ) antibiotics. Remarkably, similar SDRs are widely spread in bacteria and able to inactive other THIQ members including the clinical anticancer drug, ET-743. These findings not only fill in the missing intracellular events of temporal-spatial shielding mode for cryptic self-resistance during THIQs biosynthesis, but also exhibit a sophisticated damage-control in secondary metabolism and general immunity toward this family of antibiotics.

Versatile selective evolutionary pressure using synthetic defect in universal metabolism.

The non-natural needs of industrial applications often require new or improved enzymes. The structures and properties of enzymes are difficult to predict or design de novo. Instead, semi-rational approaches mimicking evolution entail diversification of parent enzymes followed by evaluation of isolated variants. Artificial selection pressures coupling desired enzyme properties to cell growth could overcome this key bottleneck, but are usually narrow in scope. Here we show diverse enzymes using the ubiquitous cofactors nicotinamide adenine dinucleotide (NAD) or nicotinamide adenine dinucleotide phosphate (NADP) can substitute for defective NAD regeneration, representing a very broadly-applicable artificial selection. Inactivation of Escherichia coli genes required for anaerobic NAD regeneration causes a conditional growth defect. Cells are rescued by foreign enzymes connected to the metabolic network only via NAD or NADP, but only when their substrates are supplied. Using this principle, alcohol dehydrogenase, imine reductase and nitroreductase variants with desired selectivity modifications, and a high-performing isopropanol metabolic pathway, are isolated from libraries of millions of variants in single-round experiments with typical limited information to guide design.

Dynamic behavior of liquid droplets with enzyme compartmentalization triggered by sequential glycolytic enzyme reactions.

Dynamic droplet formation via liquid-liquid phase separation (LLPS) is believed to be involved in the regulation of various biological processes. Here, a model LLPS system coupled with a sequential glycolytic enzymatic reaction was developed to reproduce the dynamic control of liquid droplets; (i) the droplets, which consist of poly-L-lysine and nucleotides, compartmentalize two different enzymes (hexokinase and glucose-6-phosphate dehydrogenase) individually, accelerating the overall reaction, and (ii) each enzymatic reaction triggers the formation, dissolution and long-term retention of the droplets by converting the scaffold nucleotides. This model system will offer a new aspect of enzymes associated with LLPS in living cells.

From a dimer to a monomer: Construction of a chimeric monomeric isocitrate dehydrogenase.

Many isocitrate dehydrogenases (IDHs) are dimeric enzymes whose catalytic sites are located at the intersubunit interface, whereas monomeric IDHs form catalytic sites with single polypeptide chains. It was proposed that monomeric IDHs were evolved from dimeric ones by partial gene duplication and fusion, but the evolutionary process had not been reproduced in laboratory. To construct a chimeric monomeric IDH from homo-dimeric one, it is necessary to reconstitute an active center by a duplicated region; to properly link the duplicated region to the rest part; and to optimize the newly formed protein surface. In this study, a chimeric monomeric IDH was successfully constructed by using homo-dimeric Escherichia coli IDH as a start point by rational design and site-saturation mutagenesis. The ~67 kDa chimeric enzyme behaved as a monomer in solution, with a Km of 61 µM and a kcat of 15 s-1 for isocitrate in the presence of NADP+ and Mn2+ . Our result demonstrated that dimeric IDHs have a potential to evolve monomeric ones. The evolution of the IDH family was also discussed.

Catalytic Control of Spiroketal Formation in Rubromycin Polyketide Biosynthesis.

The medically important bacterial aromatic polyketide natural products typically feature a planar, polycyclic core structure. An exception is found for the rubromycins, whose backbones are disrupted by a bisbenzannulated [5,6]-spiroketal pharmacophore that was recently shown to be assembled by flavin-dependent enzymes. In particular, a flavoprotein monooxygenase proved critical for the drastic oxidative rearrangement of a pentangular precursor and the installment of an intermediate [6,6]-spiroketal moiety. Here we provide structural and mechanistic insights into the control of catalysis by this spiroketal synthase, which fulfills several important functions as reductase, monooxygenase, and presumably oxidase. The enzyme hereby tightly controls the redox state of the substrate to counteract shunt product formation, while also steering the cleavage of three carbon-carbon bonds. Our work illustrates an exceptional strategy for the biosynthesis of stable chroman spiroketals.